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This  paper  addresses  the  interactions  between  homogeneous  and  heterogeneous  reactions  for  different 
hydrocarbons  namely,  compressed  natural  gas  (CNG),  liquefied  petroleum  gas  (LPG),  butane  and 
dimethyl  ether  (DME)  over  platinum.  Experiments  are  performed  to  study  the  effects  of  varying  the  tem¬ 
perature  of  the  incoming  mixture  (Tjet),  its  equivalence  ratio  (0)  and  the  Reynolds  number  (Re),  on  the 
reactivity  limits.  Computational  fluid  dynamic  (CFD)  calculations  using  detailed  chemical  kinetics  for 
both  the  platinum  surface  and  gas  phase  are  completed  for  a  range  of  methane-air  mixtures  to  resolve 
the  impact  of  varying  Tjet.  0  and  Re  on  the  compositional  structure  of  the  flow.  Comparison  between 
numerical  and  experimental  results  is  performed  where  relevant. 

For  flameless  conditions  (defined  by  the  presence  of  reactions  on  the  plate  without  a  gaseous  flame),  it 
is  found  for  all  fuels  studied  here  that  the  temperature  of  the  platinum  plate,  resulting  from  reactions 
with  the  co-flowing  fuel-air  mixture,  increases  with  increasing  7)ct  and  Re.  However,  with  CNG,  the  tem¬ 
perature  of  the  plate  peaks  near  stoichiometry  while  for  LPG,  butane  and  DME  the  peak  occurred  at  richer 
mixtures  of  0  ss  1.5.  The  reactive  limits  for  CNG,  propane  and  DME  are  found  to  broaden  significantly. 
Numerical  simulations  show  very  good  agreement  with  the  measured  plate  temperature  at  different 
equivalence  ratios.  The  computed  compositional  structure  confirms  the  existence  of  a  flame  inhibition 
effect  due  to  the  presence  of  the  catalyst  and  shows  interesting  trends  of  some  species  at  different  Re, 
Tjet  and  0.  Gas  and  surface  chemistries  seem  to  affect  a  few  species  such  as  CO,  C02,  H2,  and  H20  depend¬ 
ing  on  the  conditions  of  the  co-flowing  mixture.  Minor  species  such  as  CH3,  CH20,  0,  HCO,  and  OH  are 
largely  controlled  by  gas-phase  chemistry. 

©  201 1  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

It  is  well  known  that  the  development  of  micro-combustion  de¬ 
vices  as  sources  of  power  is  attractive  but  extremely  challenging. 
The  attraction  lies  in  the  potential  to  provide  efficient,  long-life, 
miniature  sources  of  heat  that  exploit  the  high  energy  density  of 
hydrocarbon  fuels.  The  heat  generated  could  then  be  converted 
to  (i)  electrical  energy  to  power  portable  electronic  devices  and 
hence  replace  bulky  and  short-life  batteries  [1,2],  (ii)  mechanical 
power  such  as  micro-engines  [3,4]  and  micro-turbines  [5-8]  to 
propel  micro-satellites  [9,10],  or  (iii)  chemical  energy  in  processes 
such  as  fuel  reforming  or  the  decomposition  of  ammonia  to  gener¬ 
ate  hydrogen  for  fuel  cells  [11,12].  The  challenges  in  bringing  such 
concepts  to  fruition  are  huge  and  involve  cross-disciplinary  issues 
that  span  many  areas  of  research.  The  use  of  catalysts  in  micro¬ 
reactor  devices  is  common  to  enable  better  control  of  peak  temper¬ 
atures  and  facilitate  thermal  management.  A  number  of  extensive 
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reviews  covering  various  aspects  of  the  evolution  of  micro¬ 
combustion  and  micro-power  generation  devices  have  recently 
become  available  [1,2,13,14], 

This  paper  is  concerned  with  a  specific  but  critical  aspect  of 
catalytic  micro-combustion,  namely  the  interaction  between  gas 
and  surface  reactions  and  the  reactivity  of  platinum  catalysts  with 
different  hydrocarbon  fuels.  This  hetero-homogeneous  chemistry 
interplay  has  been  studied  earlier  by  various  research  groups 
[15-22]  with  a  focus  on  four  fuels,  namely  hydrogen,  carbon  mon¬ 
oxide,  methane  and  propane.  The  studies  of  Stefandis  and  Vlachos 
[23]  were  performed  using  a  channel  flow  sandwiched  between 
two  parallel  catalytic  plates  and  focused  on  propane  (C3H8).  Other 
groups  investigated  the  behavior  of  different  fuels  such  as  hydro¬ 
gen  (H2)  [24]  and  methane  (CH4)  [25],  In  Stefanidis  and  Vlachos 
[23]  the  bulk  inlet  velocity  of  the  (lean)  fuel-air  mixture  was  varied 
as  well  as  the  separation  between  the  plates  and  their  conductivi¬ 
ties.  It  was  concluded  that  while  low  inlet  velocities  favor  surface 
reactions,  homogeneous  gas  chemistry  continued  to  play  a  role, 
albeit  decreasing,  with  decreasing  gap  size,  to  separations  as  low 
as  200  pm.  Mantzaras  and  his  group  [21]  studied  a  similar  burner 
geometry  and  fuels  but  with  the  recent  addition  of  CO/H2-air 
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mixtures.  Their  detailed  LIF-Raman  measurements  of  species  con¬ 
centrations  using  LIF-Raman  scattering  confirm  their  calculations 
and  show  that  the  reactor  stability  improves  at  higher  pressure 
[22], 

Two  other  key  findings  emerge  from  the  studies  of  Karagianni- 
dis  et  al.  [26].  The  Lewis  number  effects  are  significant  such  that 
methane  with  a  Le  number  close  to  unity  shows  “higher  transverse 
transport  towards  the  catalytic  surface”  [26]  and  hence  wider  sta¬ 
bility  limits  than  propane.  The  other  finding  is  that  the  surface 
chemistries  of  CO  and  H2  tend  to  decouple  at  high  temperatures 
while  for  lower  temperatures  around  710-720  K,  CO(s)  coverage 
of  the  surface  tends  to  inhibit  reaction  [21].  The  competition  be¬ 
tween  surface  and  gas  chemistries  near  the  catalytic  surface  is  dri¬ 
ven  to  a  large  extent  by  what  Deutschmann  et  al.  [27]  refer  to  as 
the  “sticking  probability”  of  species  on  the  surface.  The  decreasing 
sticking  probabilities  of  H2,  02,  and  CH4  explain  the  contrasting 
behavior  of  CH4-02,  and  H2-02  mixtures  on  a  platinum  surface 
as  the  percentage  of  oxygen  in  the  mixture  increases  [27]. 

It  is  evident  from  these  studies  that  in  catalytic  combustion, 
there  are  numerous  controlling  parameters  that  include  the  veloc¬ 
ity,  temperature  and  composition  of  the  gas  mixture  as  well  as  the 
properties  of  the  reactor  which  include  the  critical  dimensions 
where  reactions  occur,  surface  quality  and  thermal  properties.  This 
paper  complements  research  reported  to  date  by  presenting  a  com¬ 
prehensive  and  systematic  study  of  the  effects  relating  to  the  gas 
mixture  namely  the  effects  of  temperature,  velocity  and  equiva¬ 
lence  ratio  of  the  fuel-air  mixture  as  key  controlling  parameters 
affecting  reaction  with  the  platinum  catalyst. 

Other  aspects  of  novelty  in  the  work  presented  here  include:  (i) 
the  experimental  investigation  of  four  hydrocarbon  fuels,  namely 
compressed  natural  gas  (CNG),  liquefied  petroleum  gas  (LPG), 


commercial  butane  and  dimethyl  ether  (DME),  and  (ii)  the  use  of 
computational  fluid  dynamics  (CFD)  with  detailed  chemical  kinetic 
mechanisms  for  the  gas  and  solid  phases  to  compute  the  reactive 
structure  for  methane-platinum  and  compare  with  experimental 
measurements  where  relevant.  The  broadening  in  the  reactive  lim¬ 
its  of  these  fuels  over  a  range  of  inlet  mixture  temperatures  is  con¬ 
sistent  with  the  findings  reported  earlier  for  propane  fuels  in  the 
Swiss  roll  burner  [28],  It  should  be  noted  that  the  configuration 
studied  here  is  not  confined  but  rather  a  free  flowing  mixture  of 
fuel-air  over  a  platinum  surface.  This  is  similar  to  the  configuration 
studied  by  Kyritsis  and  co-workers  [29]  who  investigated  a 
platinum  plate  preheated  with  a  methane/air  mixture  leading  to 
a  flameless  reaction  which  was  sustained  on  the  surface.  Studies 
of  surface-gas  chemistry  interactions  in  open,  unconfined  flows 
such  as  those  reported  here  and  in  Ref.  [29]  are  still  very  relevant 
to  the  case  of  confined  micro-reactors. 

2.  Experimental  and  numerical  issues 

2.1.  Experimental 

A  schematic  of  the  experimental  apparatus  is  shown  in  Fig.  1. 
The  30  cm  long  heater  (TEMPCO)  consists  of  a  ceramic  shell  con¬ 
taining  a  heating  coil  with  a  diameter  of  8  cm  surrounding  a 
24  mm-OD  and  23  mm-ID  stainless  steel  tube.  The  fuel-air  mix¬ 
ture  flows  over  baffle  plates  positioned  inside  the  stainless  steel 
tube  to  increase  the  heat  transfer  and  achieve  exit  temperatures 
up  to  600  °C  for  Reynolds  numbers  ranging  from  250  to  1500  and 
with  equivalence  ratios  that  range  from  0  =  0  to  0  =  oo.  An  arrange¬ 
ment  of  meshes  is  used  within  the  tube  and  a  sintered  bronze  plate 
is  placed  at  the  exit  plane  to  keep  the  flow  exiting  the  tube  laminar 
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Fig.  1.  Schematic  of  the  experimental  set  up  showing  the  platinum  plate  engulfed  by  gases  from  the  jet  (i.e.  being  within  the  “valid  cone”). 
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(Re  <2000)  and  uniform  in  both  temperature  and  velocity.  Gas 
issuing  from  the  23  mm-ID  stainless  steel  tube  flows  over  both 
sides  of  the  platinum  plate,  which  has  a  width  of  6  mm,  length  of 
20  mm  and  a  thickness  of  0.25  mm.  This  region  over  which  flow 
from  the  tube  extends  without  being  affected  by  the  laboratory 
air  is  referred  to  here  as  “the  valid  region”  and  this  is  shown  sche¬ 
matically  in  Fig.  1. 

Temperature  measurements  using  a  thermocouple  were  made 
at  various  distances  from  the  exit  plane  and  above  the  platinum 
plate  to  determine  the  extent  of  the  region  over  which  the  fuel 
mixture  is  unaffected  by  the  laboratory  air.  Mapping  this  “valid  re¬ 
gion”  is  important  to  ensure  that  the  entire  platinum  plate  is  ex¬ 
posed  to  fuel-air  mixture  issuing  only  from  the  23  mm-ID 
stainless  steel  tube.  The  extent  of  the  “valid  region”  depends  on 
the  Reynolds  number,  Red  (based  on  the  burner  diameter)  of  the 
mixture  and  it  was  found  that  for  the  platinum  plate  to  be  fully 
within  this  "valid  region”  a  Reynolds  number  of  more  than  200 
should  be  used.  All  measurements  reported  here  ensured  that  the 
platinum  plate  is  fully  engulfed  by  the  mixture  issuing  from  the 
jet  so  that  the  results  are  neither  dependent  on  the  separation  from 
the  leading  edge  of  the  nozzle  nor  on  the  nozzle  diameter.  A  PID 
controller  (Novus  N1200)  is  used  to  regulate  the  temperature  of 
the  mixture  at  the  exit  plane  of  the  burner  which  is  measured 
using  a  type-K  thermocouple.  The  temperature  of  the  gases  exiting 
the  tube  can  be  varied  within  the  range  from  ambient  up  to  600  °C. 
Tylan  and  Bronkhorst  mass  flow  controllers  are  used  to  control  the 
flow  rates  of  fuels  and  air,  respectively. 

Four  different  fuels,  namely  CNG,  LPG,  commercial  butane  and 
dimethyl  ether  (DME)  have  been  tested  over  platinum  at  various 
conditions.  The  parameters  varied  during  the  experiments  are  (i) 
the  equivalence  ratio,  0  of  the  fuel  mixture  co-flowing  over  the 
platinum,  (ii)  its  temperature  at  the  jet  exit  plane,  Tjet  and,  (iii) 
its  Reynolds  number  Red  which  is  based  on  the  23  mm-ID  of  the 
fuel  tube.  Table  1  provides  details  of  the  composition  of  each  par¬ 
ent  fuel  as  obtained  from  the  supplier  and  lists  for  each  fuel  the 
range  covered  for  each  of  the  three  controlling  parameters. 

Once  steady  state  is  established  for  a  particular  flow  condition, 
and  if  the  fuel-air  mixture  does  not  self-ignite  on  the  platinum  sur¬ 
face,  ignition  is  induced  using  a  blow  torch.  If  a  gaseous  flame  is 
maintained  immediately  after  ignition  then  the  flame  is  forced  to 
blow  off.  Subsequent  to  this,  one  of  the  following  scenarios  may 
then  prevail: 

•  Extinction:  The  platinum  plate  cannot  sustain  reaction  and  its 
temperature  stabilizes  at  that  of  the  co-flowing  mixture. 

•  Flaming  combustion:  Subsequent  to  the  initial  ignition  and 
blow-off,  a  flame  may  re-ignite,  propagate  back  and  stabilize 
on  the  exit  plane  of  the  fuel  tube. 

•  Flameless  combustion:  No  flame  is  re-ignited  subsequent  to  the 
initial  blow-off  but  the  platinum  plate  remains  re-active.  While 
in  some  cases  the  plate  may  be  clearly  red-hot,  this  does  not 
have  to  be  the  case  since  surface  reactions  may  be  maintained 
at  much  lower  temperatures  as  will  be  evident  from  the  results 
presented  here. 

The  surface  temperature  of  the  platinum  was  measured  using  a 
two-color  infrared  pyrometer  (Mikron  M90R)  that  covers  a  range 


between  700  °C  and  2000  °C.  For  surface  temperatures  below 
700  °C,  an  estimate  is  obtained  by  using  a  type-K  thermocouple 
to  measure  the  gas  temperature  just  downstream  the  platinum 
plate. 

2.2.  Numerical 

Calculations  are  presented  here  only  for  methane.  The  commer¬ 
cial  package  FLUENT-6.3.26  [30]  is  used  to  perform  two- 
dimensional  calculations  for  the  experimental  conditions  studied 
here  with  the  numerical  domain  shown  in  Fig.  2.  The  domain  ex¬ 
tends  for  1 5  mm  in  the  x  direction  and  6  mm  in  the  y  direction 
and  has  been  meshed  using  8600  quadratic  grid  cells.  Monitoring 
the  species  as  well  as  the  temperature  profiles  resulting  from  a 
range  of  grids  (not  shown  here)  confirm  that  this  mesh  provides 
a  grid-independent  solution.  Special  care  is  taken  near  the  surface 
because  reactions  start  and  stabilize  there.  To  avoid  high  aspect 
ratios,  the  pave  method  is  adapted  for  meshing.  Assuming  symme¬ 
try  on  the  platinum  plate,  only  half  of  the  domain  is  calculated 
with  one  side  bounded  by  the  plate  and  the  other  by  the  fuel/air 
mixture.  The  inlet  boundary  condition  is  set  to  velocity  inlet  with 
a  flat  velocity  profile  to  match  the  experimental  conditions.  The 
plate  has  a  thickness  of  25  mm  but  only  half  of  it  (12.5  mm)  is 
modeled  because  of  symmetry.  The  leading  edge  of  the  plate  is  lo¬ 
cated  1  mm  above  the  jet  exit  plane.  The  boundary  condition  at  the 
trailing  edge  of  the  plate  (downstream)  is  modeled  as  a  pressure 
outlet.  This  is  justified  with  such  configurations  where  flow  rever¬ 
sal  may  occur  within  the  few  first  time  steps.  The  boundary  oppo¬ 
site  the  plate’s  surface  was  set  to  symmetry  since  it  is  far  away 
from  the  reaction  zone.  The  overall  extent  of  the  numerical  domain 
(6  mm  x  15  mm)  is  suitable  since  further  extensions  were  found  to 
yield  no  change  in  the  computed  flame  structure. 

The  DRM22  [31]  mechanism  is  used  for  the  volumetric  reac¬ 
tions  while  Deutschmann’s  [32]  mechanism  is  employed  to  repre¬ 
sent  surface  chemistry.  The  GRI2.11  transport  database  file  is  used 
along  with  a  comprehensive  thermodynamic  database  accounting 
for  the  relevant  gaseous  and  surface  species.  Since  the  flow  is  lam¬ 
inar,  non-unity  Lewis  number  is  adopted  with  full  multi-compo¬ 
nent  diffusion.  Gas  chemistry  is  implemented  using  the  ISAT 
algorithm  with  an  ISAT  error  tolerance  of  le-6  which  has  been 
demonstrated  earlier  to  be  quite  adequate  for  calculations  involv¬ 
ing  methane  [33].  Surface  chemistry  is  implemented  using  direct 
integration.  Radiation  heat  losses  at  the  platinum  plate  are  com¬ 
puted  using  the  Discrete  Ordinates  (DO)  radiation  model.  It  is 
known  that  the  emissivity  of  platinum  (e)  varies  with  surface  tem¬ 
perature  and  increases  from  e  =  0.1  at  a  temperature  of  550  °C  to 
£  =  0.19  at  1500  °C  [34].  A  fixed  value  for  the  emissivity,  £  =  0.15 
which  corresponds  to  a  temperature  of  1200  °C  is  adopted  here 
for  convenience.  This  is  considered  to  be  adequate  since  1200  °C 
is  close  to  the  average  temperature  experienced  by  the  platinum 
during  the  experiments  reported  here.  It  should  be  noted  that,  as 
the  emissivity  increases,  the  computed  peak  temperature  de¬ 
creases  but  the  overall  shape  of  the  computed  profiles  remains 
unchanged.  The  calculations  for  radiation  are  performed  at  every 
flow  iteration. 

The  following  computational  process  is  followed  for  all  the 
cases  presented  in  this  paper:  For  a  given  set  of  parameters 


Table  1 

Details  of  the  volumetric  fuel  compositions  and  the  experimental  conditions  covered  in  this  paper. 

Fuel  Tjet  0  Rea  Volumetric  fuel  composition 

CNG  21-400  °C  0.1-OO  250-1500  88.8%  CH4,  7.8%  C2H4,  1.9%  C02  and  1.2%  N2  with  the  remaining  0.3% 

being  a  mixture  of  propane,  propene,  butane  and  pentane 

LPG  21-400  °C  O.l-oo  250-1500  95%  C3H8, 4%  C4H10,  1%  C5  +  hydrocarbons 

Commercial  butane  21-400  °C  O.l-oo  250-1500  95.3%  C4H10, 1.9%  C3HS,  1.3%  C4HS  and  1.3%  C5H12 

DME  21-400  °C  O.l-oo  250-1500  99%  DME 
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Fig.  2.  Computational  domain  showing  dimensions  as  well  as  the  platinum  plate 
used  as  a  plane  of  symmetry. 

(equivalence  ratio,  temperature  and  Reynolds  number),  a  solution 
is  first  obtained  for  the  non-reacting  case,  where  the  flow,  energy 
and  species  transport  equations  are  solved.  This  provides  the  start¬ 
ing  flow  and  mixing  fields  to  initiate  calculation  for  the  reacting 
cases  where  the  platinum  plate  is  now  patched  with  a  1227  °C 
(1500  l<)  temperature  and  the  chemical  reactions  (volumetric  and 
surface  depending  on  the  need)  are  switched  on.  Note  that  the  en¬ 
ergy  equation  is  solved  within  the  solid  plate  so  that  heat  conduc¬ 
tion  can  be  modeled.  Second  order  discretization  schemes  are 
utilized  for  all  the  equations  and  the  under  relaxation  parameters 
are  modified  to  help  achieve  convergence  to  a  stable  solution. 
The  calculations  are  transient  and  the  time-step  used  is  initially 
1  ps  but  this  is  subsequently  increased  to  10  ps  after  the  initial 
transients  are  washed  out.  After  the  transient  solution  is  almost 
converged,  the  steady  solver  is  switched  on  for  a  few  iterations 
to  reach  a  steady-state  solution  that  is  subsequently  used  for  com¬ 
parison  with  experiments. 

3.  Experimental  results 

Measurements  of  the  peak  temperature  of  the  platinum  surface, 
Tpmax  are  presented  for  the  four  fuel  mixtures  studied  here,  namely 
CNG,  LPG,  butane  and  DME.  Figures  3  and  4  show  profiles  of  the 
peak  plate  temperature  (TPmax)  versus  equivalence  ratio  (note  the 
logarithmic  scale)  plotted  for  Reynolds  numbers  ranging  from 
250  up  to  1500.  Three  plots  are  shown  for  each  fuel  corresponding 
to  initial  approach  temperatures  of  21  °C,  200  °C  and  400  °C.  Open 
symbols  on  the  figures  refer  to  flameless  conditions  while  filled 
symbols  refer  to  flaming  combustion.  Solid  lines  refer  to  measure¬ 
ments  of  the  platinum  plate  made  using  the  pyrometer.  Dashed 
lines  refer  to  conditions  below  the  threshold  of  the  pyrometer  so 
gas  temperatures  are  measured  using  a  thermocouple  positioned 
just  above  the  platinum  plate.  Note  however  that  the  extent  of 
the  plots  shown  in  Figs.  3  and  4  do  not  necessarily  indicate  the  full 
reactive  limits.  These  are  shown  in  full  and  further  discussed  in 
Fig.  5.  The  horizontal  arrow,  where  present,  marks  self  ignition  of 
the  fuel  mixture  on  the  platinum  plate.  The  absence  of  an  arrow 
indicates  that  ignition  was  forced  using  the  torch  as  described 
earlier. 


Figure  3  shows  results  for  CNG  and  LPG  and  Fig.  4  reports  mea¬ 
surements  in  butane  and  DME.  While  there  are  some  similarities 
that  are  summarized  later  in  the  text,  each  of  the  four  fuels  shows 
sufficient  differences  to  warrant  a  separate  discussion.  The  follow¬ 
ing  observations  can  be  made  for  CNG  (LHS  plots  of  Fig.  3): 

•  Regardless  of  the  temperature  of  the  incoming  CNG-air  mixture 
Tjet,  the  maximum  value  of  peak  plate  temperature  occurs  very 
close  to  stoichiometric  and  increases  with  increasing  Reynolds 
number.  At  Tjet  =  21  °C,  the  peak  measured  platinum  tempera¬ 
ture  at  0  =  1  increases  from  Tp  m ax  =  1260  °C  to  1460  °C  as  the 
Reynolds  number  increases  from  250  to  1500.  At  higher  inlet 
temperatures  of  Tjet  =  400  °C,  the  corresponding  increase  in 
Ip ,max  at  0  =  1  is  a  bit  higher  and  ranges  from  1500  °C  at 
Re  =  250  to  1670  °C  at  Re  =  1400. 

•  The  lean  reactive  limit  decreases  with  the  increasing  tempera¬ 
ture  of  the  incoming  fuel  mixture  from  0  =  0.8  at  Tjet  =  21  °C 
to  0  =  0.5  at  Tjet  =  400  °C.  This  decrease  is  independent  of  the 
Reynolds  number  (except  for  a  slight  dependence  at 
Tjet  =  21  °C).  It  should  be  noted  that  for  gaseous  CNG-air  mix¬ 
tures,  the  lean  reactive  limits  are  0.48  and  0.35  at  mixture  tem¬ 
peratures  of  21  °C  and  400  °C,  respectively.  This  delayed 
reactivity  for  lean  mixtures  in  the  presence  of  platinum  is  also 
reproduced  numerically  (as  is  discussed  later  in  this  paper) 
and  is  due  to  the  inhibiting  effects  of  oxygen  as  it  sticks  to  the 
platinum  surface. 

•  For  rich  mixtures,  the  decrease  in  peak  plate  temperature 
occurs  in  two  stages  depending  on  the  equivalence  ratio.  A 
region  of  fast  decay  in  Tp  max  that  extends  from  stoichiometric 
to  an  intermediately  rich  equivalence  ratio  where  Tpmax 
decreases  to  around  850  °C.  This  limit  of  Tpmax  ~  850  °C  is 
reached  around  0  =  1.7  at  Tjet  =  21  °C  regardless  of  Re.  However, 
at  higher  values  of  Tjet.  a  dependence  on  Reynolds  number  is 
observed  so  that  at  Tjet  =  400  °C  the  900  °C  limits  are  reached 
at  0  =  1.7  at  Re  =  250  and  0  =  3.0  at  Re  =  1400.  For  richer  mix¬ 
tures,  the  decrease  in  Tpjnax  is  minimal  (region  of  slow  decay 
in  Tp  max)  and  seems  to  be  largely  independent  of  the  Reynolds 
number  and  even  Tjet-  This  is  expected  since  these  fuel-rich  con¬ 
ditions  are  dominated  by  surface  chemistry  and  the  surface 
sites  are  expected  to  be  saturated  with  carbon  atoms,  C(s). 

The  behavior  of  LPG  over  platinum  is  significantly  different 
from  CNG  as  can  be  clearly  seen  from  the  plots  presented  on  the 
RHS  of  Fig.  3.  The  following  observations  are  made,  noting  differ¬ 
ences  and  similarities  between  LPG  and  CNG  fuels: 

•  Unlike  CNG,  and  for  flameless  conditions,  the  tendency  for  LPG- 
air  mixtures  over  platinum  is  to  reach  peak  temperatures  at  a 
rich  (rather  than  stoichiometric)  equivalence  ratio  of  around 
0  =  1.5.  This  is  true  regardless  of  the  Reynolds  number  or  even 
the  temperature  of  the  incoming  mixture  as  is  evident  from 
the  profiles  shown  at  Tjet  =  21  °C  and  from  the  only  flameless 
condition  obtained  at  Tjet  =  200  °C  and  Re  =  250. 

•  As  marked  by  the  solid  circles  on  the  RHS  of  Fig.  3,  flaming  com¬ 
bustion  occurs  for  LPG-air  mixtures  at  Tjet  =  200  °C  and  Re  =  750 
and  1400  and  for  all  Reynolds  numbers  at  Tjet  =  400oC.  Note 
here  that,  as  expected,  flaming  combustion  occurs  within  the 
flammability  limits  of  the  gaseous  mixtures  only  and  flameless 
conditions  are  resumed  outside  these  limits.  Also,  when  flaming 
combustion  occurs,  the  peak  plate  temperature  is  dictated  by 
the  gas  temperature  and  is  independent  of  Reynolds  number 
since  the  stain  rates  are  very  low  in  this  flow  configuration.  Val¬ 
ues  of  Tpmax  =  2070  °C  and  2160  °C  are  measured  at  stoichiom¬ 
etric  for  Tjet  =  200  °C  and  400  °C,  respectively.  Note  also  that 
these  flaming  conditions  have  not  been  observed  in  the  CNG 
flames. 
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Fig.  3.  Maximum  measured  surface  temperature  plotted  versus  equivalence  ratio,  0  (on  a  logarithmic  scale)  for  CNG  and  LPG  fuels  at  different  values  of  Tjet  and  Re. 


•  As  indicated  by  the  arrow  on  some  RHS  plots  of  Fig.  3,  self¬ 
ignition  occurs  for  certain  mixtures  depending  on  the  tempera¬ 
ture,  equivalence  ratio  and  Reynolds  number.  This  is  another 
difference  with  CNG-air  mixtures  where  self-ignition  did  not 
occur  at  any  of  the  tested  conditions.  For  LPG-air  mixtures  with 
Tjet  =  400  °C  self-ignition  occurs  for  mixtures  with  0  >  0.1 
regardless  of  the  Reynolds  numbers  while  for  Tjet  =  200  °C  self 
ignition  occurs  at  0  >  0.9  regardless  of  the  Reynolds  number. 

•  A  point  of  similarity  with  CNG-air  is  that,  for  flameless  condi¬ 
tions,  the  peak  platinum  temperature  increases  with  Reynolds 
number  (but  the  peak  occurs  at  a  different  equivalence  ratio 
for  CNG  and  LPG  fuel  mixtures).  For  LPG-air  at  Tjet  =  21  °C,  val¬ 
ues  of  Tp  max  measured  at  the  peak  condition  of  0=1.5  range 
from  1310  °C  at  Re  =  250  to  1520  °C  at  Re  =  1500. 

Figure  4  shows  profiles  of  T/>max  versus  equivalence  ratio,  0  for 
butane  and  DME  plotted  at  similar  conditions  to  those  shown  in 
Fig.  3  for  CNG  and  LPG.  The  common  feature  for  all  four  fuels  is 
the  fact  that,  for  flameless  conditions,  the  peak  surface  tempera¬ 
ture,  TPmax  increases  with  Reynolds  number.  This  feature  is  consis¬ 


tent  with  earlier  studies  of  micro-reactors  performed  by  Stefanidis 
and  Vlachos  [23]  and  is  caused  by  the  increased  platinum  reactiv¬ 
ity  due  to  the  injection  of  additional  fuel  at  higher  Reynolds  num¬ 
bers.  While  this  is  one  of  the  few  common  features  noted  between 
the  four  fuels,  many  similarities  are  observed  in  the  trends  of  LPG, 
butane  and  DME  (but  not  CNG)  as  noted  here: 

•  In  flameless  combustion  modes,  peak  surface  temperatures 
occur  at  rich  equivalence  ratios  of  around  0  &  1.5  for  all  three 
fuels.  This  is  clear  evidence  that  the  platinum  surface  chemistry 
is  controlling  at  these  conditions. 

•  Flaming  combustion  occurs  at  selectively  higher  temperatures 
depending  on  the  fuel-mixture  and  the  Reynolds  numbers.  For 
butane  and  DME,  a  single  condition  of  flaming  combustion 
occurs  at  7]et  =  21  °C,  Re  =  1500  and  an  equivalence  ratio  of  1.4 
and  1.2  for  butane  and  DME,  respectively.  However,  for 
Tjet  =  200  °C  butane  sustains  a  flame  only  at  Re  =  1 500  within 
the  range  1.1  <  0  <  1.7  while  for  DME  ignites  a  flame  at  all  Rey¬ 
nolds  numbers  albeit  with  a  decreasing  range  of  lean  equiva¬ 
lence  ratios  (the  lean  limits  for  flaming  combustion  are  1.1,  1, 


822 


J.A.  Badra,  A.R.  Masri /  Combustion  and  Flame  159  (2012)  817-831 


Butane  at  21  °C  DME  at21°C 


<l>  <t> 


Fig.  4.  Maximum  measured  surface  temperature  plotted  versus  equivalence  ratio,  0  (on  a  logarithmic  scale)  for  butane  and  DME  fuels  at  different  values  of  Tjet  and  Re. 


and  0.9  at  Re  =  250,  750  and  1400,  respectively).  At  Tjet  =  400  °C, 
both  butane  and  DME  stabilize  flames  at  all  Reynolds  numbers 
and  within  a  fixed  range  of  equivalence  ratios  which  are 
0.8  <  0  <  1 .7  for  butane  and  0.7  <  0  <  2  for  DME. 

•  As  for  LPG-air,  self  ignition  occurs  at  high  inlet  temperatures  for 
both  butane  and  DME.  However,  it  is  clear  that  DME  is  much 
more  reactive  than  both  butane  and  LPG.  As  Tjet  rises  to 
200  °C,  self  ignition  of  DME-air  mixtures  occurs  at  all  Reynolds 
numbers  for  0  >  0.1.  This  is  also  true  for  Tjet  =  400  °C.  Another 
indicator  of  the  higher  reactivity  of  DME  is  that  at  rich  mixtures, 
the  transition  from  the  region  of  fast  decay  to  the  region  of  slow 
decay  in  Tp  max  moves  from  0  =  2  at  7)et  =  21  °C  to  0  =  3  at  200  °C 
and  400  °C. 

The  reactive  limits  for  the  four  fuels  studied  here  are  pre¬ 
sented  in  Fig.  5  as  plots  of  equivalence  ratio  versus  Reynolds 
number  for  various  fuel  mixture  temperatures,  Tjet  (except  for 
DME  where  only  one  temperature  of  Tjet=21°C  is  presented). 
Note  that  both  axes  show  logarithmic  scales.  Lean  and  rich  reac¬ 
tive  limits  (0L,  0R)  are  presented  for  each  value  of  Tjet.  It  should  be 


noted  that  these  limits  are  determined  here  from  temperature 
measurements  obtained  using  either  the  pyrometer  (for  the  tem¬ 
perature  range  700  °C  to  2000  °C),  or  a  thermocouple  placed  on 
top  of  the  platinum  plate.  Therefore,  reported  temperatures  out¬ 
side  the  range  of  700  °C  to  2000  °C  are  for  the  gas  phase  rather 
than  the  platinum  surface.  The  reactive  limits  are  defined  by 
the  condition  when  the  temperature  rise  decreases  to  almost 
zero,  i.e.  when  the  temperature  of  the  platinum  (or  the  gas  down¬ 
stream  of  it,  as  measured  by  the  thermocouple)  drops  to  that  of 
the  incoming  fuel  mixture  stream. 

For  all  fuels  studied  here,  the  reactive  limits  are  generally 
broader  than  those  obtained  for  the  gas-phase  alone  and  may  ex¬ 
tend  to  extreme  conditions  depending  on  the  fuel.  Selected  lean 
and  rich  reactive  limits  obtained  experimentally  for  four  fuels  are 
shown  in  Table  2  for  a  range  of  temperatures  and  for  two  Reynolds 
numbers,  namely  Re  =  250  and  1000.  Also  shown  in  Table  2  for 
comparison  are  the  corresponding  gaseous  lean  and  rich  flamma¬ 
bility  limits  (0L,gasphase  and  0R-gasPhase)-  These  are  either  obtained 
from  the  literature  or  estimated  using  a  correlation  by  Zabetakis 
[35].  The  following  additional  observations  are  noted: 
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Fig.  5.  Plots  of  0  versus  Re  on  a  log-log  scale  showing  the  experimental  lean  and  rich  reactive  limits  of  different  fuels-air  mixtures  over  platinum  and  for  different  values  of  Tjet. 


Table  2 

Experimental  lean  and  rich  flammability  limits  obtained  on  the  platinum  surface  for  the  four  fuel-air  mixtures  at  different  temperatures  and  for  Re  =  250  and  Re  =  1000.  Also 
shown  here  are  the  limits  corresponding  to  the  gaseous-only  mixtures. 


Mixture 

Temperature  (°C) 

®L, experiment 

Re  =  250 

0L,experiment 

Re  =  1000 

^L.gasphase 

0R, experiment 

Re  =  250 

0R, experiment 

Re  =  1000 

^R.gasphase 

CNG-air 

21 

0.9 

0.8 

0.48 

35 

28 

1.4 

200 

0.6 

0.6 

0.42 

250 

63 

1.6 

400 

0.5 

0.5 

0.35 

350 

100 

1.8 

LPG-air 

21 

0.4 

0.4 

0.5 

40 

8 

2.3 

200 

0.2 

0.2 

0.44 

100 

23 

2.6 

400 

0.1 

0.1 

0.37 

150 

37 

2.9 

Butane-air 

21 

0.3 

0.3 

0.56 

2.5 

1.6 

2.6 

200 

0.1 

0.1 

0.49 

3 

1.6 

2.9 

400 

0.1 

0.1 

0.41 

5 

1.7 

3.3 

DME-air 

21 

0.2 

0.3 

0.49 

CXD 

20 

3.8 

•  The  lean  reactive  limits,  0L  are  generally  independent  of  the 
Reynolds  number  and  decrease  with  increasing  temperature 
of  the  mixture.  This  statement  is  true  for  gas  phase  mixtures 
and  appears  to  imply  that  gas-phase  chemical  kinetics  play  a 
significant  role  in  determining  0L.  However,  CNG  fuel  behaves 
differently  compared  to  the  other  fuels  where,  0L,  in  the 
presence  of  platinum,  is  higher  than  0L  for  gaseous  CNG-air 
mixtures  as  shown  in  Table  2.  On  the  other  hand,  LPG,  butane 


and  DME  show  the  opposite  trend  where,  0L  dictated  by  the 
platinum  is  always  less  than  the  corresponding  gaseous  lean 
flammability  limit  (0L,gasphase)- 

•  The  rich  reactive  limits,  0R  vary  wildly  between  different  fuels 
as  also  noted  from  the  graphs  shown  in  Fig.  5  and  from  the 
values  shown  in  Table  2  for  representative  Reynolds  numbers 
of  250  and  1000.  It  is  clear  also  from  the  values  of  0R  shown 
in  Table  2  that,  except  for  butane,  these  are  much  higher  than 
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Fig.  6.  Maximum  measured  surface  temperature  plotted  versus  equivalence  ratio,  0  (on  a  logarithmic  scale)  for  methane-air  mixtures  over  platinum:  Numerical  predictions 
for  “S  +  G”  and  “S  only”  compared  with  experimental  data  for  various  flow  conditions. 


the  rich  limits  computed  for  gaseous  mixtures  of  the  same  fuel 
and  at  the  same  temperature.  These  results,  at  least  for  LPG,  are 
consistent  with  the  findings  of  Maruta  et  al.  [28]  and  Ahn  et  al. 
[36]  who  reported  a  significant  increase  in  the  reactive  limits 
for  propane  when  platinum  is  introduced  in  the  Swiss  Roll  reac¬ 
tor.  The  proximity  of  the  gaseous  and  surface  0R  limits  for 
butane  needs  further  study  and  may  be  due  to  the  platinum  sur¬ 
face  becoming  saturated  with  carbon  atoms,  C(s)  so  that  gas¬ 
eous  chemistry  dominates. 

•  The  trend  observed  for  CNG,  LPG  and  DME  (not  shown  here)  is 
for  0R  to  increase  at  higher  values  of  Tjet  regardless  of  the  Rey¬ 
nolds  number.  Conversely,  butane  shows  a  somewhat  contra¬ 
dictory  behavior,  where  0R  is  independent  of  Tjet  at  high 
Reynolds  numbers  (Re  >500).  This  is  even  inconsistent  with 
the  trends  obtained  for  the  gaseous  butane-air  mixtures  (see 
Table  2)  and  such  behavior  requires  further  investigation.  How¬ 
ever,  at  low  Reynolds  numbers  (Re  =  250)  0R  of  butane  increases 
with  Tjet. 

•  At  low  Reynolds  numbers  (say  less  than  300),  values  of  0R  start 
off  generally  a  bit  higher  and  decrease  slowly  to  plateau  at 
Re  ~  700  and  higher.  This  is  true  for  all  fuels  shown  here  except 
DME  where  the  increase  in  0R  at  low  Reynolds  numbers 
(Re  <  400)  is  somewhat  exponential  since  pure  DME  at  21  °C 


reacts  on  platinum  at  Re  =  250.  This  threshold  increases  to 
Re  =  800  when  DME  is  heated  to  200  °C  (not  shown  here).  These 
findings  support  the  view  that  DME  is  the  most  reactive  of  the 
four  fuels  tested  here. 

4.  Numerical  results 

To  confirm  the  validity  of  the  numerical  approach,  this  section 
starts  by  comparing  the  measured  and  computed  temperatures 
for  a  selection  of  cases.  Note  that  CNG  is  conveniently  used 
throughout  as  equivalent  to  pure  methane.  Given  that  boundary- 
layer  thickness  of  about  1  mm  closer  to  the  leading  edge  is  compa¬ 
rable  to  the  width  of  the  plate,  the  two-dimensionality  of  the 
boundary  layer  may  be  questioned.  To  alleviate  this  concern,  an 
experiment  was  conducted  with  a  wider  strip  of  20  mm  and  this 
reproduced  very  similar  results  to  the  6  mm  strip  within  the  region 
of  interest  reported  here  (within  at  least  5  mm  from  the  leading 
edge).  This  confirms  that  the  boundary  layer  may  be  treated  as 
two-dimensional  and  hence  the  2D-calculations  reported  here 
are  well  justified.  Next,  some  details  of  the  compositional  structure 
are  presented.  For  each  case,  calculations  are  repeated  with  the 
gaseous  reactions  being  disabled  so  that  only  platinum  surface 
chemistry  is  active.  Such  cases  are  labeled  as  “Surface”  (or  with 


J.A.  Badra,  A.R.  Masri / Combustion  and  Flame  159  (2012)  817-831 


825 


subscript  S)  as  opposed  to  the  standard  calculations  that  employ 
both  gas  and  surface  chemistry  (labeled  with  a  subscript  S  +  G,  or 
Surface  +  Gas).  Such  dual  calculations  highlight  the  importance  of 
the  homogeneous  and  heterogeneous  chemistries  in  the  meth¬ 
ane-platinum  cases  studied  here. 

Figure  6  shows  a  comparison  between  the  measured  and  com¬ 
puted  maximum  plate  temperature,  Tpmax  for  four  methane-air 


cases,  namely:  Tjet  =  27  °C  and  400  °C,  Re  =  500  and  1500.  For  each 
case,  two  calculations  are  presented  for  the  surface  (S)  as  well  as 
surface  +  gas  (S  +  G)  conditions.  Results  are  plotted  versus  equiva¬ 
lence  ratio,  0  (on  a  logarithmic  scale).  All  calculations  show  no  dif¬ 
ference  between  the  computed  temperatures  from  the  S  and  S  +  G 
cases  implying  that  surface  reactions  are  dominant  at  least  for  the 
main  heat  release.  Profiles  of  computed  species  concentration 
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Fig.  7.  Contours  of  temperature,  species  mass  fractions  of  H20,  CO,  OH  and  CH20  for  (a)  Tjet  =  27  °C,  0  =  0.9  and  Re  =  500,  (b)  Tjet  =  27  °C,  0  =  0.9  and  Re  =  1 500,  (c)  Tjet  =  400  °C, 
0  =  0.9  and  Re  =  500,  (d)  Tjet  =  400  °C,  0  =  0.9  and  Re  =  1500,  (e)  TJet  =  27  °C,  0  =  2  and  Re  =  500,  (f)  Tjet  =  27  °C,  0  =  2  and  Re  =  1 500,  (g)  Tjet  =  400  “C,  0  =  2  and  Re  =  500,  and  (h) 
Tjet  =  400  °C,  0  =  2  and  Re  =  1500,  where  both  surface  and  gaseous  reactions  are  enabled. 
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shown  later  highlight  some  differences  and  hence  some  contribu¬ 
tion  of  the  gas-phase  chemistry. 

Agreement  with  the  measured  peak  surface  temperatures  is 
also  very  good  particularly  for  methane-air  mixtures  with  equiva¬ 
lence  ratios  less  than  0  —  1.4.  The  lean  reactive  limits,  which  are 
higher  than  those  of  gaseous  methane-air  mixtures  (see  Table  2), 
are  also  well-reproduced  confirming  that  the  surface  chemistry  is 
inhibiting  reactions.  The  reasons  for  such  inhibition  at  these  lean 
conditions  are  due  to  the  oxygen’s  affinity  (or  sticking  probability 
[27])  to  the  platinum  being  stronger  than  methane.  This  saturates 
the  platinum  surface  with  oxygen  hence  leaving  no  vacant  sites  for 
the  methane  to  adsorb  and  react.  At  richer  equivalence  ratios, 
calculations  show  a  continuous  decrease  in  TPimax  and  fail  to  repro¬ 
duce  dual  regions  of  fast  and  slow  decay  in  Tp  m3X  (with  a  transition 
at  0  ~  1.7)  which  has  been  observed  experimentally.  Another  dis¬ 
crepancy  lies  in  the  extent  of  the  rich  reactive  limits  where  at 
Tjet  =  27  °C,  calculations  yield  an  extinguished  condition  at  0  —  3 
while  experiments  show  a  rich  limit  of  0  —  35  (see  Fig.  5).  This  dif¬ 
ference  may  also  be  as  significant  at  higher  values  of  T)et  although  it 


is  harder  to  see  from  the  results  shown  in  Fig.  6.  Such  discrepancy 
may  well  be  due  to  the  extent  of  surface  coverage  of  vacant  sites 
(Pt(s))  with  carbon,  C(s). 

Shown  in  Fig.  7,  in  false  colors,  are  computed  contours  of  tem¬ 
perature  and  the  mass  fractions  of  H20,  CO,  OH  and  CH20  repre¬ 
senting  the  spatial  distribution  over  almost  the  entire  region 
adjacent  to  the  platinum  plate.  Calculations  include  both  surface 
and  gas  chemistries  (S  +  G)  and  the  range  of  values  computed 
for  each  scalar  are  indicated  on  the  color  bar  on  the  far  right  of 
the  page.  Each  image  covers  a  spatial  distance  of  16  mm  (verti¬ 
cally)  and  6  mm  horizontally.  The  platinum  plate  corresponds  to 
the  vertical  axis  with  the  leading  edge  at  1  mm  above  the  base 
and  the  horizontal  axis  represents  distance  orthogonal  to  the 
plate.  Calculations  are  shown  for  two  values  of  Tjet  (27  °C  and 
400  °C),  Reynolds  number  (Re  =  500  and  1500)  and  equivalence 
ratios  (0  =  0.9  and  2.0).  Horizontal  rows  of  images  correspond 
to  a  particular  computed  parameter  (temperature  and  species 
mass  fractions)  while  columns  refer  to  a  particular  set  of  condi¬ 
tions  (7jet,  0  and  Re). 


Fig.  8.  Transverse  profiles  of  temperature  and  species  mass  fractions  computed  for  CH4-air  over  platinum  at  0  =  0.9  and  for  a  range  of  values  of  Tjet  and  Re  as  indicated  in  the 
inset. 
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A  common  feature  to  almost  all  the  images  shown  in  Fig.  7  is 
that  reaction  is  generally  initiated  within  at  the  leading  edge  of 
the  platinum  and  the  boundary  layer  formed  on  the  plate  grows 
in  thickness  with  distance.  Just  downstream  of  the  leading  edge, 
the  temperature  and  composition  fields  become  uniform  and  do 
not  change  further  with  axial  distance.  It  is  also  notable  that  the 
computed  peak  mass  fraction  of  water  (second  row  of  images)  re¬ 
mains  almost  constant  regardless  of  the  equivalence  ratio,  Re  and 
Tjet-  While  the  mass  fraction  of  CO  (third  row  of  images)  is  almost 
unaffected  by  changes  in  Re  and  Tjet,  it  increases  significantly  for 
rich  mixtures  with  peak  values  of  YCo  changing  from  0.0002  at 
0  =  0.9  to  0.146  at  0  =  2.0  for  7Jet  =  27  °C  and  Re  =  500.  The  mass 
fraction  of  OH  is  affected  somewhat  by  equivalence  ratio,  Re  and 
Tjet  but  the  changes  are  hard  to  quantify  from  the  shown  images 
and  are  discussed  later  in  this  paper.  Formaldehyde,  CH20  in¬ 
creases  in  magnitude  with  increasing  Re  and  Tjet*  For  richer  mix¬ 
tures  (0  =  2),  the  initiation  of  CH20  is  delayed  downstream  of  the 
leading  edge  of  the  plate  but  the  peak  mass  fraction  of  CH20  also 
increases  with  rising  Tjet. 


Figures  8  and  9  show  further  details  of  the  compositional  struc¬ 
ture  within  the  boundary  layer  at  two  equivalence  ratios  of  0  =  0.9 
and  0  =  2,  respectively.  Each  figure  shows  transverse  profiles  for 
temperature  and  the  mass  fraction  of  eleven  species  (CH4,  02,  CO, 
C02,  H20,  H2i  OH,  CH20,  O,  CH3,  HCO)  plotted  along  a  line  5  mm 
from  the  leading  edge  of  the  plate  (this  location  is  marked  in 
Fig.  7a).  This  position  is  deemed  to  be  representative  of  the  entire 
boundary  layer  as  can  be  seen  from  the  contours  shown  in  Fig.  7. 
The  calculations  shown  here  employ  both  surface  and  gas  chemis¬ 
tries  (S  +  G).  Each  plot  of  Figs.  8  and  9  shows  three  profiles  plotted 
versus  distance  from  the  surface  of  the  plate  and  each  profile 
corresponds  to  a  different  combination  of  Reynolds  number  and 
Tjet  as  shown  in  the  legend. 

The  temperature  profiles  confirm  earlier  findings  showing  that 
the  temperature  of  the  plate  increases  with  increasing  Tjet  and  Re 
reaching  a  value  of  1600  °C  at  Tjet  =  400  °C  and  Re  =  1500.  This  is 
well  above  the  ignition  temperature  of  methane  and  yet  no  flaming 
combustion  is  established  because  both  02  and  CH4  are  consumed 
on  the  platinum  to  produce  H20  and  C02.  This  is  evident  from  Fig.  8 


Fig.  9.  Transverse  profiles  of  temperature  and  species  mass  fractions  computed  for  CH4-air  over  platinum  at  0  =  2  and  for  a  range  of  values  of  Tjet  and  Re  as  indicated  in  the 
inset. 
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which  shows  for  0  =  0.9  mass  fractions  of  about  0.1  and  0.15  for 
H20  and  C02,  respectively.  The  rate  of  consumption  of  02  and 
CH4  increases  with  increasing  Tjet  and  Re  as  reflected  in  the  steeper 
slopes  and  this  phenomenon  of  flame  inhibition  due  to  the  pres¬ 
ence  of  platinum  is  consistent  with  what  has  already  been  reported 
by  Bui  et  al.  [37].  It  is  worth  noting  that  for  rich  mixtures  (0  =  2), 
and  at  Tjet  =  27  °C  and  Re  =  500  (Fig.  9),  the  platinum  surface  be¬ 
comes  saturated  so  that  a  residue  mass  fraction  of  methane 
1'ch4  =  0.03  is  now  seen  on  the  surface  and  the  conversion  effi¬ 
ciency  to  H20  and  C02  is  reduced  as  reflected  by  the  decrease  in 
their  respective  mass  fractions. 

The  behavior  of  other  minor  species  shown  in  Figs.  8  and  9  is 
very  interesting.  For  lean  mixtures  of  0  =  0.9  (Fig.  8),  radicals  such 
as  O,  OH,  CH3  and  HCO  form  very  close  to  or  on  the  surface  and  de¬ 
cay  very  quickly,  within  about  0.5  mm  from  the  surface  to  yield 
peaks  in  CO,  H2  as  well  as  formaldehyde  CH20.  This  behavior  is 
most  noticeable  for  Tjet  =  400oC  and  Re  =  1500  with  the  peaks 
decreasing  at  Re  =  500  and  decreasing  much  further  at  Tjet  =  27  °C 
and  Re  =  500  so  that  at  this  latter  condition  the  conversion  to 
H20  and  C02  is  at  its  highest  state.  Formaldehyde  maintains  its 


presence  on  the  surface  yet  peaks  at  some  distance  of  about 
0.5  mm  at  Ijet  =  400  °C.  This  is  driven  by  the  varying  production/ 
consumption  rates  of  CH20  which  start  negative  on  the  surface 
(consumption)  but  change  to  positive  (production)  0.5  mm  away 
from  the  surface.  Increasing  the  temperature  from  Tjet  =  27  °C  to 
400  °C  (at  Re  =  500)  changes  the  consumption  rate  at  the  surface 
by  an  order  of  magnitude,  from  -0.012  kg/m3  s  to  -0.12  kg/m3  s, 
respectively.  However,  at  0.5  mm  from  the  surface,  the  production 
rate  of  CH20  increases  by  20  times  from  0.002  kg/m3  s  to  0.04  kg/ 
m3  s  hence  explaining  the  higher  mass  fractions  of  formaldehyde 
at  this  condition. 

For  rich  mixtures  with  0  =  2  (Fig.  9),  the  abundance  of  methane 
close  to  the  surface  causes  higher  site  occupancy  and  hence  poorer 
conversion  to  stable  species  such  as  H20  and  C02.  Instead,  CO  and 
H2  are  now  much  more  abundant  than  for  the  lean  case  (0  =  0.9, 
Fig.  8)  and  peak  right  on  the  surface.  The  peak  mass  fraction  of 
C02  is  also  reduced  drastically  from  around  0.15  for  0  =  0.9  to  less 
than  0.0003  for  0  =  2.  An  interesting  shift  in  trend  is  also  observed 
for  CH20  and  HCO  in  Fig.  9,  where  the  mass  fractions  of  those  spe¬ 
cies  increase  with  increasing  Tjet  and  reduce  with  increasing  Re. 


Fig.  10.  Transverse  profiles  of  temperature  and  species  mass  fraction  computed  for  Tjet  =  27  °C,  0  =  0.9  and  Re  =  500.  Profiles  are  shown  for  two  downstream  locations  (0  mm 
and  5  mm)  using  surface  only  or  surface  +  gas  reactions  as  marked  on  the  inset. 
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The  remaining  minor  species  shown  in  Fig.  9,  (OH,  0  and  CH3) 
experience  significant  change  when  Tjet  increases  from  27  °C  to 
400  °C  but  remain  constant  with  increasing  Reynolds  numbers. 

4.1.  Homogeneous/heterogeneous  reactions 

Calculations  presented  in  Fig.  6  show  that  the  peak  surface  tem¬ 
perature  of  the  platinum  surface  Tp  max  is  mainly  controlled  by  sur¬ 
face  chemistry  since  it  remains  unchanged  when  the  gas  chemistry 
is  numerically  switched  off.  This  section  provides  a  more  thorough 
analysis  of  the  interplay  between  surface  and  gas  reaction  and 
their  impact  on  the  compositional  structure  of  the  flow  near  the 
platinum  surface.  Calculations  using  surface  only  (S)  and  sur¬ 
face  +  gas  (S  +  G)  chemistry  are  repeated  for  lean  and  rich  methane 
flows  with  0  =  0.9  and  0  =  2,  at  two  inlet  temperatures,  Tjet  =  27  °C 
and  400  °C  and  two  Reynolds  numbers  (Re  =  500  and  1500). 

Figures  10-12  show  sample  results  for  only  three  cases  with 
Re  =  500  and  (i)  0  =  0.9,  rjet  =  27°C,  (ii)  0  =  0.9,  rjet  =  400°C,  and 
(iii)  0  =  2,  Tjet  =  27  °C.  Each  figure  shows  profiles  of  the  computed 
temperature  and  mass  fraction  of  eleven  species  (CH4,  02,  CO, 
C02,  H20,  H2,  OH,  CH20,  O,  CH3,  HCO)  plotted  transverse  to  the 


platinum  plate  starting  from  two  locations,  the  leading  edge 
(0  mm)  and  5  mm  from  leading  edge.  The  locations  of  these  pro¬ 
files  are  shown  by  black  arrows  on  the  top  LHS  of  Fig.  7.  They  are 
chosen  here  to  illustrate  the  spatial  evolution  of  the  boundary  layer 
noting  that  the  location  of  5  mm  from  the  leading  edge  of  the  plate 
is  representative  of  the  remaining  boundary  layer  (see  Fig.  6).  For 
each  species,  and  for  each  of  the  two  locations  on  the  plate  two 
profiles  are  shown  for  the  (S  +  G)  and  (S)  calculations. 

It  is  shown  from  Fig.  10(0  =  0.9,  Tjet  =  27  °C,  and  Re  =  500)  that 
temperature  and  the  computed  mass  fractions  of  CH4,  02,  CO,  C02 
and  H20  are  almost  fully  controlled  by  surface  chemistry  given  that 
the  profiles  from  S  and  S  +  G  calculations  almost  overlap  (except 
for  a  slight  deviation  in  the  CO  profiles).  However,  H2,  OH,  CH20, 
O,  CH3  and  HCO  show  huge  differences  confirming  that  gas  chem¬ 
istry  plays  a  key  role  in  their  formation.  This  is  expected  for  CH20, 
O  and  HCO  which  are  not  part  of  the  surface  mechanism  anyway. 
The  methyl  radical,  CH3  is  part  of  the  platinum  chemistry  but  re¬ 
acts  only  on  the  surface  and  never  desorbs  as  a  gaseous  species. 
Hydrogen  mass  fraction  is  zero  on  the  plate’s  surface  and  peaks 
some  0.5-1.0  mm  away  but  is  formed  by  virtue  of  homogeneous 
chemistry.  Conversely,  OH  mass  fraction  at  the  leading  edge  is 


Fig.  11.  Transverse  profiles  of  temperature  and  species  mass  fraction  computed  for  Tjet  =  400  °C,  0  =  0.9  and  Re  =  500.  Profiles  are  shown  for  two  downstream  locations  (0  mm 
and  5  mm)  using  surface  only  or  surface  +  gas  reactions  as  marked  on  the  inset. 
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Fig.  12.  Transverse  profiles  of  temperature  and  species  mass  fraction  computed  for  Tjet  =  400  °C,  0  =  2  and  Re  =  500.  Profiles  are  shown  for  two  downstream  locations  (0  mm 
and  5  mm)  using  surface  only  or  surface  +  gas  reactions  as  marked  on  the  inset. 


highest  on  the  plate’s  surface  and  this  is  almost  entirely  due  to  sur¬ 
face  chemistry  but  further  away  from  the  plate,  gas  chemistry  con¬ 
tributes  to  the  destruction  of  OH  as  is  evident  from  the  profiles 
shown  in  Fig.  10. 

This  is  confirmed  from  a  comparison  of  the  rates  of  production 
of  OH  across  the  boundary  layer.  At  5  mm  from  the  leading  edge, 
the  net  reaction  rate  for  OH  changes  from  1.1  e~8  kg/m3  s  to  0  kg/ 
m3  s  across  a  2.5  mm  boundary  layer  using  surface  chemistry 
alone.  With  surface  +  gas  chemistries  at  the  same  location  the  net 
reaction  rates  increases  form  -0.018  kg/m3  s  to  0kg/m3s  across 
a  0.25  mm  boundary  layer.  This  confirms  that  the  OH  which  des¬ 
orbs  from  the  platinum  by  a  set  of  surface  reactions  gets  consumed 
by  the  gaseous  reactions  as  it  diffuses  away  from  the  surface  into 
the  boundary  layer,  it  should  be  noted  that  the  pattern  described 
in  Fig.  10  for  the  case  of  0  =  0.9,  Tjet  =  27  °C,  and  Re  =  500  is  also  rel¬ 
evant  at  higher  Reynolds  numbers  with  0  =  0.9,  7)et  =  27  °C,  and 
Re  =  1500  (not  shown  here). 

Figure  1 1  shows  computed  profiles  for  the  case  where  Tjet  rises 
to  400  °C  at  0  =  0.9  and  Re  =  500.  The  overall  structure  is  very  sim¬ 
ilar  to  the  case  shown  in  Fig.  10  for  Tjet  =  27  °C  except  for  CO  and 


C02  where  the  gas  chemistry  has  more  impact.  The  computed  mass 
fractions  of  CO  have  increased  significantly  for  the  S  +  G  case  at 
5  mm  from  the  leading  edge  of  the  platinum.  The  higher  tempera¬ 
ture  of  the  incoming  mixture  activates  CO  production  due  to  gas¬ 
eous  reactions  some  0.5  mm  away  from  the  plate  where  CH20  is 
also  formed  in  sizeable  quantities.  The  computed  profiles  of  C02 
profiles  show  also  a  slight  deviation  due  to  the  active  gas 
chemistry. 

For  a  richer  incoming  mixture  with  0  =  2,  Tjet  =  400  °C,  and 
Re  =  500,  Fig.  12  shows  significant  changes  in  the  compositional 
structure  compared  to  the  case  presented  in  Fig.  11  for  0  =  0.9. 
The  mass  fractions  of  CO  and  H2  formed  on  the  platinum  increase 
by  more  than  three  orders  of  magnitude.  Additionally,  the  overlap 
between  the  S  and  S  +  G  profiles  of  H2  and  CO  implies  that  surface 
chemistry  is  in  control.  The  mass  fraction  of  C02  decreases  by  a 
corresponding  amount  but  the  same  cannot  be  said  about  H20 
which  drops  only  slightly.  The  C02  levels  decrease  significantly 
and  some  active  gas  chemistry  is  evident  at  5  mm  from  the  leading 
edge  increasing  the  peak  C02  mass  fraction  from  about  0.001  to 
0.0022.  The  fact  that  H20  does  not  decrease  by  an  amount 
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corresponding  to  the  increase  in  H2  is  largely  due  to  lower  mass 
fraction  of  OH  obtained  for  this  case.  The  increase  in  the  mass  frac¬ 
tion  of  H2  is  an  artifact  of  the  increased  formation  of  H(s)  on  the 
surface  due  to  the  higher  concentration  of  CH4  in  the  mixture 
which  then  recombines  to  form  H2  through  a  series  of  surface 
reactions. 

Albeit  limited  to  methane  for  now,  it  is  clear  that  a  wealth  of 
information  can  be  gleaned  from  these  detailed  and  computation¬ 
ally  expensive  calculations.  While  gaseous  mechanisms  for  the 
remaining  fuels  (propane,  butane  and  DME)  may  be  made  avail¬ 
able,  their  detailed  surface  chemistry  on  platinum  is  severely  lack¬ 
ing.  Earlier  calculations  by  Stefandis  and  Vlachos  [23]  and 
Karagiannidis  et  al.  [26]  used  a  single-step  catalytic  reaction  for 
propane  on  platinum.  There  is  an  urgent  need,  therefore,  to  devel¬ 
op  detailed  catalytic  mechanisms  for  the  range  of  fuels  studied 
here  as  well  as  others. 


5.  Conclusions 

The  experimental  and  numerical  studies  reported  here  confirm 
that  the  hetero-homogeneous  interactions  in  the  presence  of  plat¬ 
inum  are  complex  and  highly  dependent  on  the  fuel  used  as  well  as 
other  parameters  such  as  temperature,  equivalence  ratio  and  Rey¬ 
nolds  number  of  the  gaseous  fuel-air  mixture.  The  following  con¬ 
clusions  are  drawn: 

•  The  peak  temperature  of  the  platinum  surface,  Tp  max  increases 
with  increasing  Reynolds  number  for  all  four  fuels  studied  here. 
However,  the  surface  temperature  peaks  near  stoichiometric  for 
CNG-air  mixtures  whereas,  for  LPG,  butane  and  DME  fuels, 
TPi max  occurs  at  an  equivalence  ratio  of  about  1.5. 

•  The  reactivity  limits  of  the  fuel-air  mixtures  are  generally 
broader  in  the  presence  of  platinum  except  for  butane  where 
the  lean  and  rich  reactive  limits  remained  close  to  those  of 
the  gas  phase.  DME  is  the  most  reactive,  decomposing  on  plat¬ 
inum  as  a  pure  fuel,  followed  by  CNG  and  then  LPG. 

•  Very  good  agreement  is  obtained  between  the  calculated  and 
measured  temperatures  over  a  range  of  conditions  for  CNG- 
air  mixtures  where  7]et,  Re  and  0  are  varied.  It  is  also  found  that 
the  peak  platinum  temperature  is  largely  unaffected  by  the  gas- 
phase  reactions. 

•  Combined  gas  and  surface  chemistries  affect  a  few  species  such 
as  CO,  C02,  H2,  and  H20  depending  on  the  conditions  (Tjet,  Re 
and  0)  of  the  co-flowing  mixture.  Minor  species  such  as  CH3, 
CH20,  0,  HCO,  and  OH  are  mainly  controlled  by  gas-phase 
chemistry. 
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